A SOLIDIFICATION PHENOMENON IN RANDOM PACKINGS 

L. BOWEN, R. LYONS, C. RADIN, AND P. WINKLER 

Abstract. We prove that uniformly random packings of copies of a certain simply-connected 
figure in the plane exhibit global connectedness at all sufficiently high densities, but not at 
low densities. 

1. Introduction 

The densest way to cover a large area with non-overlapping unit disks is as in Figure [TJ 




Figure 1. The densest packing of unit disks in the plane 

A packing is a collection of congruent copies of a subset with pairwise disjoint interiors. See 
|Fej| for a proof that the above packing is indeed densest possible for unit disks. 

It is an old unsolved problem to understand whether densest packings of spheres, simplices 
or other shapes, in a Euclidean or hyperbolic space of any dimension, exhibit crystallographic 
symmetry. For instance, this is the spirit of Hilbert's eighteenth problem; see |Fej[ IRadj for 
background. 

Using physics models of two- and three-dimensional matter as a guide, we are tempted 
to try to gain insight about densest packings by considering packings at densities below 
the maximum. (For an example concerning spheres in M 3 , see |KRS .) In effect, we are 
emphasizing not so much that densest packings and sparse packings differ by their symmetry, 
as that they differ in some fundamental geometric fashion. Indeed, it is commonly suggested 
in the physics literature (see for instance AHJ) that 2-dimensional models of matter do 
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not exhibit crystallographic symmetry, and it is sometimes said by mathematicians that in 
high dimensional Euclidean space, densest packings of spheres may not have crystallographic 
symmetry. So perhaps it is appropriate to reexamine the precise manner in which densest 
packings differ fundamentally from sparse packings, and to use packings at less than optimum 
density guide. 

The density of the unit disk packing of Figure [T] is d := n/ \^12 rs 0.91; packings of density 
0.89 (say) can be obtained by shrinking the disks slightly, shaking the whole collection, and 
then expanding the whole picture to recover unit-size disks. One might ask whether by 
doing this one preserves long-range correlations. Experiments and discrete models indicate 
that such long-range correlations are indeed preserved, but this has never been proved. See 
|BLRWj for background. 

In this work we prove the conjectured behavior, not for disks but for deformed disks, copies 
of a "zipper" tile designed specifically for the purpose; see Figure 121 This tile can cover the 
plane completely, in which case the packing has density 1, and is completely connected in any 
sense. What we show is that even at somewhat lower densities, the uniform random packing 
still has rich structure; in particular it has a form of connectedness associated with site 
percolation [Gri . What this means for packing large but finite boxes (with torus boundary 
conditions) is that the necessary gross irregularities of most packings at such high densities 
occur, but not in a way to disconnect the packings. Although we define "uniform random 
packing" of the plane by limits of measures on packings of finite boxes, the key to our proof 
is to examine isometry-invariant probability measures on packings of the whole plane and to 
show that the ones that maximize "degrees of freedom per tile" are unique for high densities. 

In physics one is interested in change of behavior as density changes. We show that 
at high density there is a nonzero probability of an infinite linked component, and that 
this probability is zero at low density. Thus, there are different "phases" of the packings 
[BLRWJ. (This is the same phase transition seen in continuum percolation, where one looks 
at overlapping disks with random independent centers, but our methods are quite different. 
Indeed, no such result is known for packings of disks.) 

Although we believe such a result also holds for packings of disks or of spheres — pairs 
of which would be called "linked" if sufficiently close — we are able to prove the result only 
for our tiles, which are shaped to allow three well-defined levels of pairwise separation. It 
is generally understood that crystalline behavior is not seen in two dimensions, so the form 
of connectedness we use may be useful in understanding the role of geometry in Hilbert's 
problem. 

2. Description of the tile 

We consider packings by a deformed disk denoted by t, referred to as "the tile" and depicted 
in Figure 121 In this section, we define it precisely. 

Let H be a regular hexagon of area 1. Let r be the radius of the in-circle of H. Let D 
be a disk concentric with the in-circle and of radius r + p, where < p <C 1 is a number we 
shall choose more precisely later. We shall construct the shape t by modifying H as follows; 
D will be called the shadow disk of t. 



A SOLIDIFICATION PHENOMENON IN RANDOM PACKINGS 



3 



Figure 2. The zipper tile 




Figure 3. Close-up view of the fringe 

As shown in Figure El the tile t equals H with each side modified by a "fringe" and each 
corner modified by a hook and inlet, where a hook is about half an element of the fringe. 
As shown in Figure El the fringe height is 2 p. The elements of the fringe have two different 
size "necks", one of size p 2 and one of size 2p 2 , allowing neighboring tiles to be linked in 
either of two well-defined modes, "tight linked" and "loose linked" , the former illustrated in 
Figure 0] and the latter illustrated in Figure We say that two tiles t are linked (tightly or 
loosely) if when one is held fixed, the other can be moved continuously only by a bounded 
amount (without overlapping the first). A tight link is one that permits no movement of 
one tile while fixing the other, while a link that is not tight is called loose. A key feature 
of our model is that when two tiles are tightly linked, any motion of one would necessitate a 
corresponding motion of the other. As we shall explain, the uniform probability distribution 
on packings of the plane at given density is a limit of such a distribution on packings of larger 
and larger tori. In our model, these distributions on packings of finite tori are concentrated 
on packings with the maximal number of degrees of freedom, and therefore intuitively the 
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fewest possible number of tiles bound by tight links. This gives us useful control on the 
packings in the support of our distributions. 




Figure 4. Tight-linked tiles 

A tile is called fully linked on one side if it is linked with another tile on that side in 
such a way that either they are tight linked and the line joining their centers goes through 
the midpoint of the sides of the corresponding hexagons; or the tiles are not tight linked 
but can be moved continuously so that their shadow disks touch each other. A tile is fully 
linked if it is fully linked on all sides. We note that the fully tight-linked packing (Figure HJ 
corresponds to a tiling by the original hexagon and has density 1, and that the tile has area 
1 by construction. 

3. Statement of Results 

To state our results we need some notation. Let X be the space of all packings of the 
plane by the tile t. Given a compact subset K of the plane and two packings of the plane, 
we consider the distance between the two packings with respect to K to be the Hausdorff 
distance between the unions of the tiles in the respective packings intersected with K. Then 
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Figure 5. Loose-linked tiles 



X is endowed with the topology of Hausdorff convergence on compact subsets; X is compact. 
Intuitively, two packings are close in X if they are close in the Hausdorff sense in a large ball 
centered at the origin. We shall define a probability measure on X that is "uniform" on the 
set of all packings of a fixed density. For this, we shall need the space X n of all packings by 
the tile of the n x n torus R 2 /(nZ) 2 . 

For any integer m, let X n m C X n consist of those packings which contain exactly m tiles 
(X njTn is empty if m is large enough). To each tile, we assign the set of 6 unit vectors based 
on its center and pointing to the center of each of its edges. Through this assignment, we 
can view X n ^ m as a subset of T™/H m , where T n is the unit tangent bundle of the n x n torus 
modulo a 2tt/6 rotation and the symmetric group S m acts by permuting the factors. 

When m/n 2 is small, X n>m is (3m)-dimensional. However, when m/n 2 is sufficiently large, 
the dimension of X„ im inside T™ /E m is less than 3m. This is because at least two tiles 
in any packing of X n)Tn will have to be tightly linked, so that it is impossible to move one 
continuously without moving the other. Thus it is useful to decompose X n ^ m into a (finite) 
disjoint union of sets X n ^ m ^ of packings containing exactly k tight links. Generically, the 
dimension of X n ^ m ^ is 3 (m — k). The dimension can be strictly less than this if the packings 
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are jammed in the sense of [DTSCJ, although this fact will not be important for us. The 
top dimension of X n>m means the maximum dimension of all X n)Tn ^. Let p, n>m be the 
probability measure on X n>m obtained by normalizing the Hausdorff measure on X n m in 
the top dimension of X n ^ m with respect to the natural metric inherited from T™ /H m . We 
interpret p, n ,m as being a uniform measure. The fact that p, n>m is supported on those packings 
with the fewest number of tight links will be crucial in the analysis to follow. 

Let X n be the space of all (nxn)-periodic packings of the plane. In other words, X n consists 
of those packings that are preserved under translations by nL x nL. Under the quotient map, 
this space is naturally identified with X n . Therefore, we can view the measures \i n , m as living 
on X n C X. 

For a fixed density d G [0,1], let pS d > be any measure obtained as the weak* limit of 
measures of the form p n ^ m such that n — > oo and m/n 2 — > d. (Note that m/n 2 is the density 
of every packing in the support of p, n>m and d is the average density of a packing chosen with 
respect to pf^\ see Lemma 15.11 ) A priori, ji^ may not be unique, although we shall prove 
that it is for large enough d. 

As we have seen, two tiles t can be linked. Therefore it makes sense to speak of a linked 
component of a packing; it is a maximal sub-packing such that for every two tiles t, t' in it, 
there is a sequence t = ti, t 2 , • • • , t n = t' such that ti is linked to U + i (i — 1, . . . , n — 1). A 
tight-linked component is defined similarly, but we require £, to be tightly- linked to ti + \. 

We say that a measure on the space X of packings is invariant if it preserved under the full 
isometry group of the plane. All the measures we consider are probability measures unless 
stated otherwise. 

Let Ao be the unique invariant measure on tilings (packings that cover R 2 ) by our tile. 
Let Ai be the unique invariant measure on packings by t such that all tiles are fully loose 
linked, are as close as possible to each other, and the packing has hexagonal symmetry. Write 
A s := sAi + (1 - s)X Q . 

Our main results are the following: 

Theorem 3.1. There exists < dx < 1 such that if d > d\, p^ is unique and equals \ s , 
where s := (1 — d)/(l — d\). 

Corollary 3.2. The fj,^ -probability that the origin is inside a tile belonging to an infinite 
linked component is nonzero for d > di . 

Proposition 3.3. For some c?2 > 0, the probability (with respect to any /j,^ for any d < d 2 ) 
that the origin is inside an infinite linked component is zero. 

4. Tile properties 

Lemma 4.1. For small p, if tiles t\ and ti are not tightly linked and do not overlap, then 
the distance between their centers is at least 2r + 2p. 

Proof. Consider the line segment from the center of t\ to the center of ti. If this segment 
traverses near a corner of t\ or t 2 , then it must be longer than 2r + 2p for small enough 
p. Supposes it crosses a fringe of t\ and of t 2 - If the tiles are not linked, then the claim is 
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obvious. If they are linked, then to minimize the distance, it must be that their fringes match 
up (so they are fully linked on one side). Thus the closest they can come is if the two are 
pushed flat up against each other so that their shadow disks touch. In this case, the distance 
between the centers is exactly 2r + p. □ 

We shall say that two tiles are densely loose linked if they are loose linked and their 
shadow disks touch. There is a unique invariant measure on maximally dense packings by 
congruent disks BHRS 3 . Hence the probability measure Ai that we defined earlier is the 
unique invariant measure on packings by t such that all tiles are fully and densely loose 
linked. Let d\ be the density of such a packing. 

Given a tile t in a packing P, we denote by V(t) the Voronoi cell of the center of t with 
respect to the centers of the other tiles; that is, V(t) is the open set of points closer to the 
center of t than to the center of any other tile. We denote the area of a region A of the plane 
by \A\. 

Lemma 4.2. The following holds for small enough p > 0. For any packing P, if t G P is 
a tile that has no tight links, then the area ofV(t) is least \/d\. Moreover, equality holds iff 
the configuration of tiles determining V(t) is congruent to a corresponding configuration of a 
packing in the support of \\. 

Proof. For a tile t, let H(t) denote the hexagon from which t is created. For x > 0, let H x (t) 
denote the homothetic copy 1 -^H{t) about the center of H{t). 

Suppose t is a tile of P without any tight links. Consider the rays Ri, . . . , R§ from the 
center of the hexagon H(t) through each of its 6 vertices. These rays divide the plane into 6 
sectors, Si, . . . , S e . 

By construction, if t and ti are loose linked, then \H p (t) fl H p (ti)\ = 0(p 2 ): The hexagon 
interiors do not intersect if they are parallel, while if they are not parallel, they can intersect 
only very slightly at a corner. The openings at which a corner can enter have area 0(p 2 ) as 
/'- -('. ' 

Thus, we have proved that whenever t is loose linked in the sector Si, then \V(t) fl 5j| > 
\H p {t)\/6 - Si, with <5x = 0(p 2 ) as p -> 0. 

Similarly, if t and t 2 are not linked at all, then \H 2p (t) fl H 2p (t 2 )\ = 0(p 2 ): Again, their 
interiors do not intersect if they are parallel, while if they are not parallel, they can intersect 
only very slightly at a corner. So there exists 5 2 > such that whenever t is not linked in 
the sector Si, we have \V(t) fl S{\ > \H 2p (t)\/6 — 5 2 , with 5 2 = 0(p 2 ) as p — » 0. 

Therefore, if t has no tight links but is not fully linked, then 

for some j with < j < 5. Given that Si, S 2 are of order p 2 while \H 2p (t) \ — \H p (t) \ is of order 
p, for p small enough we may conclude that > \H p {t)\ in this case. 

On the other hand, the geometry of a tile is such that for small p, if t\ and t 2 are two tiles 
loose linked to t, then ti cannot be tight linked to t 2 . Now suppose that t is fully loose linked. 
Then the Voronoi cell of the center of t is determined by six tiles t\, . . . ,t 6 all loose linked 
to t and all with the property that their shadow disks D, Di, . . . , D$ do not overlap (by the 
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previous lemma). It follows |Fej2| that |V(t)| > \H p (t)\, with equality iff each of the disks 
Di, . . . , D e touches D. But there is only one way in which this can occur (up to isometry). 
So V(t) = H p (t) in this case. This implies that the configuration t,t±, . . . ,te is congruent to 
a corresponding configuration of a packing in the support of Ai. □ 

It is easy to see that given p > 0, there exists e > such that for any finite component c of 
tight-linked tiles in any packing, the union V c of the Voronoi cells of the centers of the tiles 
of c has area at least j c + ePer c . Here j c is the number of tiles in c and Per c is the perimeter 
of the union of hexagons corresponding to c. Let e be the largest such constant. Let 5 > 
be such that the area of the Voronoi cell in the fully densely loose-linked packing equals 
1 + ePeri + 5, where Peri is the perimeter of the hexagon of a single tile. Since e = p + 0(p 2 ) 
and S = 0(p 2 ), we have: 

Lemma 4.3. For sufficiently small p, there are e, 5 > such that for any finite tight-linked 
component c, 

4 



and 



1 + Per^ + 5 
\V C \ > j c + ePer c , 

5 < e/100. 

5. High Density 

Recall that X is the compact space of all packings of the plane by the tile (with the topology 
of Hausdorff- metric convergence on compact subsets). Let M be the space of isometry- 
invariant Borel probability measures on X. For any p e M, we denote by \p\ := p(A ) the 
density of p, where A is the set of all packings P G X, one of whose tiles contains the 
origin. Since a tile is the closure of its interior, A Q is a closed set. 

Lemma 5.1. If pi G M converges to p in the weak* topology, then \pi\ converges to \p\. 

Proof. Let P denote the union of tiles in a packing, P. For any invariant probability measure 
v and any z G M 2 , we have 

H = / W } MP) = J i {MeP} MP)- 

Integrating over z in a unit-area disk, D, with respect to Lebesgue measure and using Fubini's 
theorem gives the identity \u\ = J \P(lD\ du(P). Since the function P i— > |PDD| is continuous 
on X, the lemma follows. □ 



Recall that Ao is the unique invariant measure on tilings by our tile, so that |Ao| = 1. 
Recalling that d\ is the density of a fully densely loose-linked tiling, fix a density d with 
d\ < d < 1. Let pn he the uniform measure on configuration of tiles at density d^ in an 
N x N torus, where d^ — > d as N — > oo. To prove theorem 13. 11 we shall show that the weak* 
limit of hn exists and equals X s , where s := (1 — d)/(l — di). 
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We shall use several lemmas that depend on the following notation. Given a packing 
P G X, let 

• tp be the tile of P such that the origin belongs to V(tp) (this exists as long as the 
origin is not on the boundary of a Voronoi cell), 

• Kp be the tight-linked component containing tp, 

• jp be the number of tiles in Kp, and 

• f(P) '■= 3/ 'jp if jp is finite and tp contains the origin, and otherwise. 

Thus f{P), in a sense, measures the number of degrees of freedom per tile near the origin. 

Lemma 5.2. If v is any measure in M, then J f dv(f) < 3|z/| ; with equality iff t P has no 
tight links for v -almost every packing P. 

The proof is immediate. 

Lemma 5.3. // a sequence (u n ) C M converges to v in the weak* topology, then J f dv n 
converges to f fdu. 

Lemma 15.31 is proven in a manner similar to Lemma 15.11 

Given a finite tight-linked component c, let the congruence class of c be C and let Xc C X 
be the space of all packings P for which tp exists and Kp is in C. Let X" be the space of 
all packings P with density 1, where "density" refers to the usual concept of the limit of the 
proportion of the area of P inside a large disk centered at the origin as the radius tends to 
infinity. Let X' C X be the space of all packings P such that Kp is infinite and either the 
density of P is less than 1, the density is not defined, or tp does not exist. Thus, X is the 
disjoint union of X' , X" and the collection of Xc for all C . 

Let v be any invariant probability measure with density d. Let vc be v conditioned on 
Xc, v' be v conditioned on X' and v" be v conditioned on X" . Since Ao is the only invariant 
probability measure with support in X" , we have v" = A . Thus, 

v = v {X'y + z/(X")A + < X c>c ■ 



c 



Define the density \u\ := uj(A ) as before, but for any (invariant or non-invariant) proba- 
bility measure w onX. We have 



d=\v\ = u{X') 1 1/' | + u(X") + "( X c) Wc\ 

and 



c 



I fdu = Y J »(Xc) I fdv c = J2"(Xc)-Wc\, 
J c J c 3c 

where jc is the number of tiles in C. 

Lemma 5.4. Let v 6 M and C be a finite- component class. Suppose that < s < 1 is such 
that | A s | = \vc\- Then J f d\ s > f f dvc- Moreover, equality holds only if jc = 1- 
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Proof. As in the proof of Lemma f5. 11 we have that \vq\ = J jc/\V(tp)\ dvc(P). 

First suppose that j c = 1. Then \vc\ = J l/\V(tp)\duc(P) < d\ by Lemma l4~21 This 
means that s — 1 and J f dvc = 3 \vc\ — 3 |A S | = / / d\ s . 

Now assume that j = jc > 1 and put p := Perc- By definition, 

/ d\ s — s If dXi + (1 — s) / / gL\ — s f d\i = 3sdi . 



Since Uc(f) = 3\v c \/jc — 3|A s |/jc = 3(sdi + 1 — s)/j, it suffices to show that 

^■ki-.) 

J 

which is equivalent to 

s(jdi - di + 1) > 1 . 

Now sdi + (1 — s) — \vc\ < tz - j where £ is from Lemma [4.31 Solving for s gives 



l 



j 



s > 3+£p 



whence it is enough to show that 

1 - 



(jd 1 -d 1 + l) T -^>l. 
This boils down to 

di(pe + l) > 1. 

Now, j > 1 implies that p > (7/6)Peri, where Peri is the perimeter of a single tile. Since 
e/100 > S (by Lemma f4.3|) . this implies that pe + 1 > 1 + ePeri + 5 = l/d±, proving the last 
inequality. □ 

Lemma 5.5. We have J f dv < j f d\ s for all v e M with \u\ = \X S \. Equality holds only if 

• v(Xc) = for every component class C with jc > 1 and 

• whenever v(Xc) > and jc = 1, we /love = d\. 

Proof. Recall that 

y fdv = Y,"(Xc) J fdv c . 

c 

For each component class C, let sc be defined as follows: 

• if there exists s G [0, 1] such that \uc\ = sdi + (1 — s), then set sc '■= s; 

• otherwise, set sc '■= 1. 

Let c^c := Sc-Ai + (1 — sc)X and 

a := + u(X")) X + ^ z/^cVc . 

c 
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From the previous lemma, if \vq\ > d±, then J f due < / f du>c, with equality only if jc — 1- 
If \v c \ < di, then sc — 1 and 

/ tZi/ c = < 3di = / / du c ■ 

Jc J 



Summing up, we obtain 

J fda = J>PM j fduc 
c 

> J>(X C ) j f dv C = j f 

c 



dv . 



Moreover, equality holds only if v(Xc) = for every component C with jc > 1 and \ vc\ = d\ 
whenever j c = 1. Since \uc\ > Wc\, we have 



\a\ = u(X')+u(X")+^u(X c )\u c \ 

c 

> u(X') \u'\ + v{X") + "( x c) Wc\ 



c 

= \v\ 
= |A,|. 

Since a and X s are both convex combinations of Ao and Ai, this implies that j f da < j f d\ s 
with equality iff a = X s . Thus, f f du < j f d\ s . In the equality case we must have 
j f dv = j f da = j f dX s and a = X s . This implies that v{Xq) = if jc > 1 and \vc\ = d\ 
ifj c = l. ' ' ' □ 

Lemma 5.6. Let v G M. If \v\ = \X S \, then j f du < j f dX s . Equality holds iff v — X s . 

(Informally, X s uniquely maximizes the number of degrees of freedom per tile for invariant 
measures of a fixed density.) 

Proof. The previous lemma implies j f du < j f dX s . Assume j f du = j f dX s ; then 

u = u{X')u' + u{X")X + u{X c )u c 

where C is the component of size 1 and \uc\ = d\. This gives J f du = u(Xc)3dx = J f dX s = 
3sdi. Hence u(Xc) = s. Since u' has density strictly less than 1 = |A | but \u\ = \X S \, we 
must have u(X') = 0. That is, 

u = u(X")X + u(X c )u c . 

Since u and Ao are isometry invariant, vc must also be isometry invariant. By Lemma 14.21 
Ai is the unique isometry- invariant measure with support in Xc and with density d\. Hence 
vc — Ai. This implies u = X s and the proof is finished. □ 



12 



L. BOWEN, R. LYONS, C. RADIN, AND P. WINKLER 



Proof of Theorem VJ . 1\ It is easy to see that one can pack the N x N torus in such a way that 
there is a large region of tight-linked tiles and a large region of densely loose-linked tiles, in 
such a way that the interface between the two regions has a density which approaches zero as 
N tends to infinity, and the density of the packing Pn tends to d. Let ujn be the invariant 
measure supported on isometric copies of P/v (a pull-back of P/v to the plane). Then u N 
tends to A s in the weak* topology. By Lemma 1531 this implies that f f du>N — > f f d\ s . 

Now /xjv, the uniform measure of density d^ on the N x N torus, satisfies f f dfiN > 
J f du N . This is because \i n is by definition supported on packings with the maximal number 
of degrees of freedom for the given density d^. Hence liminfTv / / d^N > liminfAr J f du^ = 
JfdXs. 

Therefore, if /Xqq is any weak* subsequential limit of (^n)n, then f f d^ > f f d\ s . But 
g?at — > d, so | /loo | = I A s I by Lemma f5. II The previous lemma now implies that fXoo = A s . □ 

Returning to the discussion of the introduction, we note that from simulations of hard 
disks, one would expect the corollary to hold even for a range of densities below di, but we 
do not know how to prove this. 

Remark on Higher Dimensions 

The basic features of our argument can be generalized to dimension 3 or higher, except for 
our use in Lemma f4.2l of |Fej2| on the minimal Voronoi region in disk packings in the plane. 
It would be of interest if this part of our proof could be replaced by an argument insensitive 
to dimension. 



6. Low Density 

In this final section, we confirm the intuition that at low densities, there will be no infinite 
loose-linked component. It is obvious that there is no infinite tight-linked component at 
densities smaller than d\. 

We begin with a lemma that holds for any tile shape (in fact, for any collection of shapes 
and sizes, as long as each can be fit into a disk of some fixed radius s, and "density" is 
interpreted as number of tiles per unit area). 

Lemma 6.1. For small enough density d, if a packing P is drawn from fx^ d \ then the proba- 
bility that the disk Br of radius R about the origin contains more than 9R 2 d tile centers goes 
to zero as R — > oo. 

Proof. Let s be the radius of the smallest disk containing the tile (in our case, s is about 
2 i/2 . 3 -3/4 . (i + 2p)) and choose 

.05 

< d< To — 2 ; 

for our zipper tiles with small enough p, d < .003 suffices. Let T be the set of tiles whose 
centers fall in Br, k := \nR 2 d~\ and £ > 9R 2 d. Letting /x^(-) denote the probability of an 
event with respect to the measure fi^ d \ we shall show that 

^(\T\=i) < 
n(d)(\T\ = k) - 7 
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for some constant 7 < 1. It then follows that 

00 l(9-TT)R 2 d\ 

fi {d) {\T\ >9R 2 d) < fi (d) {\T\ = k) V 7 ^ fc < fi (d) {\T\ = k)^ >0 

^ — ' 1 —7 

as R — > 00, as desired. 

The measure /i^ is the limit of uniform distributions of configurations on the N x N torus 
T^, in turn obtainable by choosing a sequence of n — [N 2 d\ points from the Lebesgue dis- 
tribution A on T^r as the centers of the tiles, orienting each tile independently and uniformly 
at random, and finally conditioning on no overlap. We denote by A(|T| — j) the a priori 
probability that exactly j points fall inside B R (which we take to be some fixed disk in the 
torus). 

Let $ be the event that there is no overlap among the tiles whose centers lie in Br, and 
^ the event that there is no overlap involving any tile whose center falls outside Br. Then 



^ d \\T 


\=£) A(|T 


\=£) A($| 


\T 


\=i) A(tf | 


|T 


\=e A $) 


H {d) (\T\ 


= k) A(|T| 


= k) A($ 


\T\ 


= k) A(tf | 


\T\ 


= k A $) 



and our job is to bound the three fractions on the right. 

For the first, we note that \T\ is binomially distributed in the measure A, hence 



K\T\=£) = (") (jg) l 1 -^) < (n-fc)!/(n-l)! 

A(|T|=A;) (3(#) fc (i-#r fc " £!/A;! V 1 -^ 

(t/e)* ft (n — kp k \ i J 

for large R. 

The next fraction is easy: since we may throw the first k centers into B R , then the remaining 
£— k, we have 



A($ | 


|T 




A($ | 


|T| 


= *) 



is the probability that the additional £— k centers do not cause a collision, which is at most 
1. 

For the (inverse of) the third fraction, we throw n — i centers into the region outside Br, 
then the remaining t — k. A new point, if it lands at distance greater than 2s from any 
previous point or from the disk Br, causes no new overlap; and at each stage there are fewer 
than n — k points already placed. Hence 

, 'N 2 -ir(R + 2s) 2 - (n-k)4its 2 ^ e ~ k 



A(* | 


\T\ 


= k A $) 


A(tf | 


\T 


\—£ A $) 



N 2 - nR 2 
AttsR + Atts 2 \ £ ~ k . <n o.n 



/ 9 , 47rs J R + 4vrs 2 \ , 9 A , 



for iV > i?. 
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Putting the inequalities together, we have 

^ d) (\T\=i) ( .95 y- k = e _ k 

fj,W(\T\ = k) ~ VI - I3s 2 d) 1 
where 7 := .95/(1 — 13s 2 d) < 1 by choice of d. □ 

Proposition 6.2. For some di > 0, the -probability that the origin is inside an infinite 
connected component of loosely-linked tiles is zero for d < di- 

Proof. Let d G (0, .003) be a density to be chosen later. Let P be a packing drawn from 

we aim to show that the probability that the origin is connected by a loose-linked chain of 

tiles of P to some point at distance R approaches zero as R — > 00. 

We again choose some large radius R and let T be the set of tiles of P whose centers fall 
inside the disk B R . 

Fix the positions of the tiles of P \ T (the black tiles of Figure EJ) and consider the space 
of packings having these tiles plus n tiles whose centers fall in Br. We think of this space 
as being a subset of Ti(B R ) n /S n , where Ti(B R ) is the unit tangent bundle of B R (modulo a 
27r/6 rotation to take into account the symmetries of the tile) and the symmetric group acts 
by permuting the factors. 

If ot n is the volume (in Ti(£>ij) n /£ n -space) of this space and m < n, then by packing n—m 
tiles into B R and then the remaining m in the left-over space, we have 

«« > 7 „ x a n - m — W{R - 2s) 2 -mr(2s) 2 ] m , 
( ) ml L J 

\n—m/ 

where s is, as before, the radius of the circle circumscribing a tile. This takes into account 
possible intrusion of tiles in P — T into B R , and the fact that a tile center at point x can 
exclude nearby centers but only within distance 2s of x. 

Let (3 denote the "wiggle room" of a tile t loose linked to a stationary tile t', that is, the 
3-dimensional volume of the space of positions of t; then (3 = 0(p 3 ) (but we use only that (3 
is bounded by a constant). If a packing "percolates", that is, contains a chain of loose-linked 
tiles connecting the center to the boundary of B R , let t\, . . . , t m be a shortest such chain (the 
dark grey tiles of Figure Note that m > Rj (2s). For each % > 2, the tile ti is linked to one 
of the three sides of t^_i farthest from the side of t^_i linked to tj_ 2 , and has wiggle room at 
most (3 with respect to tj_i. Accounting for the orientation of t% and allowing the remaining 
n—m tile centers to fall anywhere in B R , we have that the 3n-dimensional volume of the set 
of percolating packings is bounded by (n/3) ■ 6 ■ 3 m ~ 2 . (3 m ~ 1 ■ a n - m < 3 m ' , (3 m ~^ a n - m . 

Comparing with the lower bound for a n , we find that given \T\ = n < 9dR 2 , the probability 
of percolation is less than 



3 m f3 m - 1 a n - m 3 m (3 m - 1 n\/(n-m)\ / 27f3dR 2 

c7 n - [vr( J R-2s) 2 -n7r(2s) 2 ] m < W [(R - 2s) 2 - 3QdR 2 s 2 ] 




which goes to zero as R (thus also m) increases, for suitably chosen d. Since we know from 
Lemma f6. II that /x^(|T| < 9dR 2 ) approaches 1 as R — > 00, the proposition follows. □ 
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Figure 6. An unlikely configuration of tiles in and around Br 

A more careful argument would prove Proposition ^ . 2l for any density below 1/ (47r(2/3y / 3)) = 
.2067 + for sufficiently small p, but clearly the probability of percolation will remain for much 
higher densities than that. 

7. A Conjecture 

We have shown that high-density random packings of zipper tiles in the plane contain an 
infinite loose-linked component with positive probability, while low-density random packings 
do not. What happens in the case of ordinary disks, where there is no apparent linking 
mechanism? We believe, but cannot prove, the following 

Conjecture. Suppose fj,^ is defined as above for geometric disks of radius 1. Join two 
centers by an edge if their distance is at most 2 + e for some fixed e < 1. Then for sufficiently 
high density d below the maximum, the graph resulting from a configuration drawn from (i^ 
will contain an infinite connected component a.s. 

This statement can be shown by a standard Peierls-type argument for large e; this may be 
known already. In general, there is some parameter set of (e,d) C (0, oo) x (0, 1) for which 
there is an infinite component. For small d or for large e, the problem is quite similar to 
continuum percolation, where one connects by an edge two points of a Poisson point process 
if their distance is at most r. Because of homotheties, one may fix the intensity of the point 
process to be 1. Then there is a phase transition in r. Our problem is quite different in really 
having two parameters, but our conjecture is that there is a phase transition in d for every 
e nevertheless. 
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